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Recent PHENIX measurements of the elliptic (^2) and hexadecapole (^4) Fourier 
flow coefficients for charged hadrons as a function of transverse momentum (pr), 

o 

CN| ■ collision centrality and particle species are presented and compared with results 

from the PHOBOS and STAR collaborations respectively. The status of extensions 

H- » ■ 

\q [ to future PHENIX measurements at lower beam energies is also discussed. 
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I. INTRODUCTION 

The discovery of large elliptic flow (v 2 ) for all particle species studied at the Relativistic 
Heavy Ion Collider (RHIC), is a key piece of evidence for the creation of hot and dense 
^ . partonic matter in ultra relativistic nucleus- nucleus collisions jl, 2]. This is well supported 
by the observed agreement between differential flow measurements and calculations that 
model an essentially locally equilibrated quark gluon plasma (QGP) having a small value 
for the specific shear viscosity (r]/s) [3H6|]. After ten years of experiments at RHIC, the 
most extensive set of flow measurements for various hadrons with different masses, charges, 
^ . quark content and hadronic cross-sections became available for the first time in the history 
of heavy-ion collisions. They show that, for a given centrality, elliptic flow for all observed 
hadrons [at RHIC] scale to a single curve when plotted as Vijn q versus KE^jn^ where n q is 
the number of constituent quarks in a given hadron species and KEj- is the transverse kinetic 
energy for these hadrons 3, |8j . Such scaling is illustrated in Fig. [D where a compilation of 



elliptic flow results for identified hadrons, measured by STAR [9 12] and PHENIX QQfl 
for minimum-bias Au+Au collisions at ^/snn — 200 GeV (left panel), 10-40% midcentral 
Au+Au collisions at ^snn = 62.4 GeV (middle panel) and 0-50% central Cu+Cu collisions 
at a/J/vtv = 200 GeV (right panel). The observation that all of these data show scaling 
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indicates that the agreement between the STAR and PHENIX v 2 data is better than ~ 15%, 



otherwise, the scaling would be broken. More importantly, the observed scaling suggests that 
the bulk of the elliptic flow at RHIC energies is partonic, rather than hadronic 71. IsL Il8]. 
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FIG. 1: V2/n q as a function of KET/n q for identified particle species obtained in minimum-bias 
Au+Au collisions at sJsnn = 200 GeV (left panel), in 10 - 40% midcentral Au+Au collisions at 
yfsNN = 62.4 GeV (middle panel) and in - 50% central Cu+Cu collisions at y / J/\Fiv = 200 GeV 
(right panel). 

Lowering the collision energy and studying the energy dependence of anisotropic flow 
allows a search for the onset of the transition to a phase with partonic degrees of freedom 
at an early stage of the collision (], [l?]] . In the present work, we review the most recent flow 
measurements performed by the PHENIX collaboration at top RHIC energies, and give the 
most comprehensive comparison of these results with measurements from both PHOBOS 
and STAR collaborations. We also discuss the progress of our ongoing flow analysis at lower 
collision energies at RHIC. 



II. RESULTS 



Future progress in the extraction of the transport properties of hot and dense partonic 
matter from flow measurements at RHIC, depends strongly on further developments of 
theoretical models, as well as a better understanding of the systematic errors associated with 
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flow measurements. First comparisons of differential elliptic flow data from the top RHIC 
energy (^snn = 200 GeV) with viscous relativistic hydrodynamic simulations, demonstrate 



that a 20% uncertainty in the measured Vn 
value of specific shear viscosity ( rj/s ) [3 



eads to ~ 60-70% unsertainty in the extracted 



161 ] . The role of fluctuations and so-called 



"nonflow" correlations can be important for such measurements [20j. The next lesson is 

that the simultaneous measurement of v-i and higher harmonics such as hexadecapole flow 

(uj) and triangular flow (t> 3 ), can help to better constrain r]/s j^lJ-Q]. 

PHENIX has addressed many of these issues via a new set of measurements of charged 
□ 

hadron u 2 and [19J. These measurements were performed in the two PHENIX central arms 
(M < 0.35) relative to event planes obtained from four separate reaction-plane detectors in 
the range 1.0 < \r]\ < 3.9, see left panel of Fig. 2. Multiple event planes allow a search for 
possible A?7-dependent non-flow contributions that would influence the magnitude of t>2,4, 
which is crucial for reliable extraction of transport coefficients. 
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FIG. 2: Left panel: The azimuthal angle 4> an d pseudo-rapidity 77 acceptance of PHENIX detectors 
used for event plane reconstruction, together with the central arm acceptances for charged hadron 
measurements. Right panel: comparison of V2 vs. -/Vp art and V4 vs. N part for charged hadrons 
obtained with several reaction plane detectors for the pt selections indicated [19]. Ratios for the 
Pt range 1.2-1.6 GeV/c are shown in (c) and (d). 



The right panel of Fig. 2 compares the double differential flow coefficients V2,a(pt, N part ) 
for event-plane detectors spanning the range 1.0 < |?7| < 3.9. Within systematic errors, they 
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agree to better than ~ 5% (10%) for v 2 (v^) in mid-central collisions and approximately 
10% (20%) in central and peripheral events. This agreement indicates that the present 
measurements of v-i and v± are not affected by Arj- and pr-dependent non-flow contributions 
for pt ^ 3 GeV/c and centrality < 60% 

The recently installed time-of-flight detector (TOFw) with intrinsic timing resolution of ~ 
75 ps significantly improved the particle identification capabilities, as well as the acceptance 
of PHENIX for flow measurements of identified hadrons. Time-of-flight measurements in 
conjunction with measured momentum and flight-path length allow proton/kaon separation 
up to pt — 4.5 GeV/c via mass-squared calculations. With the combined measurements of 
mass-squared and photon yield from the Aerogel Cherenkov counter (ACCV one can extend 
the t>2 measurements of charged pions and protons up to px — 6 GeV/c [24]. This is shown 
in the left panel of Fig. [3] where measured v 2 values for charged pions , kaons and protons 
are plotted in the scaled variables v-ijuq versus KEx/n q . These scaled data indicate that 
the universal KET/n q scaling, evident at low KEt, is broken after REx/rig ^ lGeV/c 2 , 
indicating a possible change in the physics. 
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FIG. 3: V2/n q as a function of KEx/n q for identified charged hadrons obtianed in 0-60% central 
Au+Au collisions at y / J/v7v = 200 GeV [ left panel ] and v n /n q ^ 2 as a function of KET/n q for the 
same particle species [ right panel ]. 

Statistically significant measurements of V4 for identified hadrons also enable a study of 
a related scaling relation between u 4 for different particle species, as well as the ratio v±jv\. 
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The right panel of Fig. [3] shows the result of such a study in which both v 2 and V4 are plotted 
using the generalized scaling variables, v n /n^ 2 vs. KE T /n q . Further detailed studies of this 
generalized scaling for different flow harmonics, as well as the KET/n q range where sca ling 
breaks down, should help to better understand the transition from soft to hard physics [16 1 
and the process of hadronization at RHIC. 

An important step toward reducing the systematic uncertainty associated with flow mea- 
surements at RHIC, is a detailed comparison of differential flow results obtained by different 
collaborations - PHOBOS, PHENIX and STAR. One way to show such comparisons is to 
plot the respective ratios of these measurements as a function of pt and centrality, while 
ensuring that both centrality and mean pt are the same. For reference, we have used 
PHENIX flow measurements of charged hadrons with respect to the MPC (3.1 < \rj\ < 3.7) 
and RXN combined (1.0 < \rj\ < 2.8) event planes 19( from recent data obtained in 2007. 
The ratios t> 2 (PHOBOS) /t> 2 (PHENIX) are plotted in Fig H]as a function of px for several 
bins in collision centrality as indicated. The values for PHOBOS ^(pt) were obtained for 
charged hadrons with < 77 < 1.5 with respect to the second order event plane measured at 
2.05 < r] < 3.2 [25]. 
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FIG. 4: pt dependence of the ratio of PHOBOS v<i values [25(] obtained by event plane method 
v 2 {EP} (2.05 < \r}\ < 3.2) to the PHENIX v 2 measured using RXN (1.0 < \r]\ < 2.8) [left panel ] 
and MPC (3.1 < \r]\ < 3.7) [right panel ]. 



The ratio plots indicate that the difference between PHENIX/PHOBOS results are 
< 5% and do not show a significant pt dependence. The STAR i> 2 for charged hadrons from 
Au+Au collisions at ^s^n = 200 GeV measured with respect to the event plane from the 
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FTPC detector (2.5 < \rj\ < 4.0) are not available. Therefore our comparison was made 
with STAR results obtained via the standard event plane method using the central TPC 
(0 < \r)\ < 1) t> 2 {EP}, and a modified event plane method where particles within |A^| < 0.5 
around the highest vt particle were excluded for the determination of the modified event 

nn 

plane t^j-E^}; this procedure reduces some of the non- flow effects at high p T [3, HlJ. The 
ratios v 2 (STAR EP)/t> 2 (PHENIX) are plotted in the top panel of Fig [5] for 30-40% and 
40-60% midcentral Au+Au collisions at A /J/v r /v = 200 GeV. 
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FIG. 5: pt dependence of the ratio of STAR i; 2 values for charged hadrons obtained by standard 
[upper panel ] and modified [ lower panel ] event plane method w 2 {EP} using central TPC (|r/| < 1.0) 
to the PHENIX tj 2 values obtained by event plane method using MPC/RXN event planes for 30- 
40% and 40-60% midcentral Au+Au collisions. 



The difference between STAR/PHENIX v 2 values is less than 2-5% below pt ~2.5 GeV/c. 
At higher transverse momentum, STAR v 2 is systematically larger than the PHENIX v 2 and 
the ratio tends to grow with p-p reaching the value of 20% at ~5.5 GeV/c. The difference 
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in v 2 values at higher can be attributed to non-flow effects due to di-jets which are mostly 
suppressed by the rapidity gaps in the case of the PHENIX measurements. The lower panels 
of Figure [5] show the same ratios but in this case the STAR v% data were obtained by the 
modified event plane method i; 2 {EP 2 }. Here, the difference between STAR/PHENIX ratios 
at high pt, are much smaller, but still persist on the level of 5-10% at p^ — 5.5 GeV/c for 
mid-central collisions. The comparison for other centralities can be found elsewhere [26]. 
Figure [6] show the comparison of the same PHENIX t> 2 data set with STAR t> 2 results 
obtained using multi-particle methods. That is, the four particle cumulant method (left 
panel) and the Lee- Yang Zero method (right panel). 
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FIG. 6: pt dependence of the ratio of STAR vi values for charged hadrons obtained by four-particle 
cumulant method f2{4} [left panel] and Lee- Yang Zero method [right panel] using charged tracks 
from central TPC (|r/| < 1.0) to the PHENIX v<i values obtained by event plane method using 
MPC/RXN event planes. 

The ratio plots indicate that STAR t> 2 results obtained using multi-particle methods 
are smaller than PHENIX t> 2 {EP} event plane results by 10-12%. Note however that this 
difference does not show a pr dependence and can be attributed to the expected difference 
in the eccentricity fluctuations for the two- and multi-particle measurements. 
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III. COLLISION ENERGY DEPENDENCE OF AZIMUTHAL ANISOTROPY 

A large number of elliptic flow measurements have been performed by many experimental 
groups at SIS, AGS, SPS and RHIC energies over the last twenty years. However, the fact 
that these data were not obtained under the same experimental conditions, do not allow 
a detailed and meaningful comparison in most cases - the situation at RHIC is of course 
somewhat better. Experimental differences include: (a) different centrality selection, (b) dif- 
ferent transverse momentum acceptance, (c) different particle species, (d) different rapidity 
coverage and (e) different methods for flow analysis. The results from PHENIX and STAR 
indicate that the magnitudes and trends of the differential elliptic flow , t>2( Pt, centrality ), 
change very little over the collision energy range ^/snn = 62 - 200 GeV, indicating saturation 
of the excitation function for t> 2 at these energies [hJ see left panel of Fig. IIHI 
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FIG. 7: Left panel: differential t>2 vs. 1 /saT/v for charged hadrons in nucleus-nucleus collisions 
for centrality cut of 13-26% and two different mean px values. From [27J. Right panel: centrality 
dependence of event plane resolution factors for V2 ( closed symbols ) and ( closed symbols ) 
measurements using PHENIX RXN detector for Au+Au collisions at ^snn = 39 - 200 GeV. 

The figure also shows that the differential elliptic flow for charged hadrons increase by 
almost 50% from the top SPS energy of 17.3 GeV to ^s NN = 62 GeV at RHIC. This con- 
clusion is based on the comparison of PHENIX results with published results from CERES 



collaboration 28]. However, the comparison of STAR results for vi^Vt ) for charged pions 
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and protons at ^snn = 62 GeV with results from NA49 Collaboration at 17.3 GeV [29], 
leads to the conclusion that the differential v% change only by 10-15% from top SPS energy 
to RHIC lOj. This may indicate that the existing flow results at the SPS are prone to 
large systematic uncertainties which are not yet well understood. Given the fact that the 
energy density increases by approximately 30% over the range ^/snn = 62.4 - 200 GeV, the 
apparent saturation of differential t>2 at RHIC could be an indication of a softening of the 
equation of state due to the crossover transition. 

In Run 2010 the PHENIX Collaboration collected ~ 5 x 10 8 minimum-bias Au+Au 
events at tJsnn — 62.4 GeV and ~ 2 x 10 8 minimum-bias Au+Au events at ^snn — 39 
GeV. The analysis of these data is progressing very well. The right panel of Fig. II III shows 
the centrality dependence of the event plane resolution factors for t> 2 and V4 measurements 
obtained for the second order event planes from the RXN detector for the collision energies 
39, 62.4 and 200 GeV. They indicate that reliable extraction of the flow harmonics should 
be straightforward. Note that this large data set in conjunction with an improved event 
plane resolution, gives an equivalent of ~ 30-fold increase in statistics over the previous 
measurement of v% at 62 GeV [271 ] . We expect that they will allow detailed measurements 
of the differential V2 and v 4 as a function of pt and centrality for several particle species at 
both energies. 
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